Finite-size scaling and deconfinement transition in gauge theories 
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A new method is proposed for determining the critical indices of the deconfinement transition in gauge theories, 
based on the finite-size scaling analysis of simple lattice operators, such as the plaquette. A precise determination 
of the critical index v, in agreement with the prediction of the Svetitsky-Yaffe conjecture, is obtained for SU(3) 
gauge theory in (2+l)-dimension. Preliminary results for SU(2) in (3-|-l)-dimension are also given. 



1. INTRODUCTION 

Given a (d + l)-dimensional pure gauge the- 
ory undergoing a continuous deconfinement tran- 
sition at the critical temperature Tc, the order 
parameter of this transition is the Polyakov loop. 
The d-dimensional effective model obtained by in- 
tegrating out all degrees of freedom except the 
order parameter is globally invariant under the 
center of the gauge group. This effective model 
possesses only short-range interactions and there- 
fore the phase transition, when continuous, is ac- 
companied by long range fluctuations only in the 
order parameter 

According to the Svetitsky-Yaffe conjecture [Q, 
the [d -\- l)-dimensional gauge theory and the d- 
dimensional effective model, if it also displays a 
second-order phase transition, belong to the same 
universality class. Therefore, their critical prop- 
erties such as critical indices, finite-size scaling, 
correlation functions at criticality, are predicted 
to coincide. The validity of the Svetitsky-Yaffe 
conjecture has been confirmed in several Monte 
Carlo analyses [||. 

In the present work we have considered 
SU(3) in (2-hl)-dimension and SU(2) in (3-hl)- 
dimension, which belong to the universality 
classes of the 3-state Potts model in 2-dimension 
and of the Z(2) (Ising) model in 3-dimension, re- 



spectively. 

2. FINITE-SIZE BEHAVIOR AT CRITI- 
CALITY 

For a d-dimensional statistical model, a way to 
extract critical indices is to study the finite-size 
behavior of suitable operators. If a mapping be- 
tween operators in the {d -I- l)-dimensional gauge 
theory and operators in the d-dimensional effec- 
tive model is established, it is possible to exploit 
finite-size effects also in the gauge theory to deter- 
mine its critical indices. This would provide with 
a new method for the computation of the criti- 
cal indices of a pure gauge theory with a second 
order phase transition. 

The Svetitsky-Yaffe conjecture intrinsically es- 
tablishes a correspondence between the Polyakov 
loop of the gauge theory and the order parameter 
of the effective model. The mapping of the pla- 
quette operator of the gauge theory with a linear 
combination of the identity and the energy oper- 
ators of the effective model has been shown in ||] 
for SU(2) in (2-|-l)-dimension, which generates 
as effective model the 2-dimensional Ising model. 
In this case, the critical behavior of the effective 
model is known exactly, thanks to the methods 
of conformal field theory. This mapping between 
operators has allowed to study, by means of uni- 
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versality arguments, several non-perturbative fea- 
tures of gauge theories Q . 

In the present work the method is apphed to 
the computation of the critical index v of the cor- 
relation length. In the c?-dimensional statistical 
model, the finite-size behavior of the (lattice) en- 
ergy operator E is given by 



(1) 



where L is the lattice size, understood to be large 
enough, and fc is a non- universal constant. In the 
(d l)-dimensional pure gauge theory, any oper- 
ator O, invariant under gauge and center symme- 
try is expected to scale according to 



= ci I 



irrclevants 



(2) 



where / is the identity operator and e the (scaling) 
energy operator of the statistical model. Simple 
operators which satisfy the mentioned symmetry 
requirements are Wilson loops and operators like 
P{x)P^ (y), where P{x) is the Polyakov loop at 
the site x and x, y represent different sites of the 
d-dimensional space. This ansatz was introduced 
and tested in and used in 

3. NUMERICAL RESULTS 

SU(3) in (2+l)-dimension|] 

The lattice operators considered were the elec- 
tric and magnetic plaquette and P{x)P^ (y), with 
X and y neighbor sites in the 2-dimensional spatial 
lattice. All these observables have the required 
symmetry properties and can be computed with 
high accuracy in Monte Carlo simulations. Simu- 
lations were performed on lattices with temporal 
size Nt = 2 and spatial sizes — Ny = L rang- 
ing from 7 to 30; (3 was set at its critical value 
f3^{Nt = 2) = 8.155, taken from Ref. §. The 
simulation algorithm was a mixture of one sweep 
of a 10-hit Metropolis and four sweeps of over- 
relaxation consisting in two updates of (random) 
SU(2) subgroups. For each simulation 400K equi- 
librium configurations were collected, separated 
each other by 10 updating steps. The error anal- 
ysis was done by the jackknife method applied to 



-"^The results presented in this subsection have been pub- 
lished in Ref. js). 



data bins at different levels of blocking. The table 
of numerical results can be found in Ref. M . 
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Figure 1. SU(3) in (2-|-l)-dimension: electric and 
magnetic plaquette vs L^/'^~'^, where v comes 
from the multibranched fit. Similar figure has 
been obtained for the Polyakov loop correlator. 



A single multibranched fit was performed of the 
three data sets, taking electric (magnetic) plaque- 
ttes from lattices with even (odd) L, in order to 
avoid cross correlations. Polyakov loop correla- 
tor data were measured separately and are there- 
fore not correlated with the plaquette measure- 
ments. The result of the fit is i/ = 0.827(22), 
with Xrcd — ^-^^ (s^® Fig.||), to be compared with 
the exact value in the 2-dimensional 3-state Potts 
model, ly = 5/6 — 0.833 A previous Monte 
Carlo determination by the method in SU(3) 
in (2+l)-dimcnsion g gave lyMC = 0.90(20). 
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SU(2) in (3+l)-dimensiorf| 

The lattice operators considered are the elec- 
tric and magnetic plaquette. Monte Carlo simula- 
tions were performed on lattices with Nt = 2 and 
Nx — Ny — L ranging from 5 to 22, at the criti- 
cal value of /3, f3c{Nt 2) = 1.8735, taken from 
Ref. . The simulation algorithm adopted is the 
over-relaxed heat-bath js) and the error analysis 
was done as in the previous case. For each simu- 
lation we have collected so far a number of con- 
figurations ranging from 50K to BOOK, according 
to the lattice size. 
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Figure 2. SU(2) in (3-|-l)-dimension: electric 
and magnetic plaquette vs L^/'^"'^, where v comes 
from the multibranched fit. 



(odd) L has given v = 0.632(8), with x^d = 0-30 
(see Fig. §. This result is to be compared 
with the determination of v in the 3-dimensional 
Ising model. In this model, the high-temperature 
expansion method gave the following very ac- 
curate result §: = 0.63002(23). Compare 
also with the Monte Carlo determination by the 
method in SU(2) in (3+l)-dimension |0,|lO|: 
VMC = 0.630(9). 
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The multibranched fit on the electric (mag- 
netic) plaquettes taken from lattices with even 

■^The following results are preliminary and have been ob- 
tained in collaboration with C. Vena of the Dip. Fisica, 
Univ. Calabria. 



